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Abstract. We study the charmonium 2S states xp' and r)'^, and the bottomonium 2S states T' and 77^, using 
the light-front quark model and the 2S state wave function of harmonic oscillator as the approximation of 
the 2S quarkonium wave function. The decay constants, transition form factors and masses of these mesons 
are calculated and compared with experimental data. Predictions of quantities such as Br(i/;' 7?7c) are 
made. The 2S wave function may help us learn more about the structure of these heavy quarkonia. 
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1 Introduction 

Charmonium physics has long been an interesting issue as 
it is related with both the perturbative and non-perturbative 
QCD. Charmonia not only provide us with the oppor- 
tunity to investigate the interactions between the con- 
stituent quarks and the structure of quarkonia, but also 
the chance to learn and understand the QCD dynamics 
better. As exited states of charmonia, ip' ^.nd ry^ have been 
studied by many authors. The decay widths of "0' — ^ e'^e" 
and T]'^ — )■ 27 were calculated with both relativistic and 
QCD radiative corrections |lj and the result -0' — >■ e+e" 
is in agreement with experimental data. The nonrelativis- 
tic potential model [5] and the Godfrey and Isgur (GI) 
model [21i4j have achieved much success, but their predic- 
tions of the decay widths of 0' — ?> 777c {iv'c) ^I'fi larger 
than experimental data. The lattice QCD result [51[3] of 
J li) ^ jT]c is consistent with experimental data, but the 
result of t/j' — >■ 777^ has too large uncertainties. The in- 
termediate meson loop contribution to the decays tp' ~^ 
IVc ilVc) was investigated recently [311], and the results 
are closer to experimental data. We also investigate these 
decays, using light-front formalism and the harmonic oscil- 
lator wave functions as the approximate wave functions of 
the IS and 2S quarkonia. In fact, there are still some puz- 
zles concerning tp' , such as the well-known "p7r puzzle" [HI 
[TDlfTTlfT^ . and the recent unanticipated small experimen- 
tal value of Br(7/>' jT])/Br{ip' 777') ^MM- ^ot the 
"pTT puzzle" , Ref. [T3] suggested the explanation that the 
T/^'pTT coupling is suppressed due to the mismatch between 
the nodeless wave function of the cc in the \udcc) Fock 
state of p and the one-node 2S cc wave function of ip' , 
and our postulation of the 2S wave function of 7/;' may be 



email: mabq@pku.edu.cn 



able to offer a numerical realization for this explanation. 
The EES and CLEO collaborations have conducted many 
experimental measurements on IS and 2S charmonia, and 
the decay mode ip' — 777^ is being studied by BES-IIL It 
is then important to learn carefully about the structures 
and decay mechanisms of the 2S charmonia. 

The 2S bottomonia have been studied by some exper- 
iments, but many data about them are still not available, 
such as the mass and decay data of 77^, [16]. In Ref. [1], 
the decay modes rj'i^ 2j and T' — ^ e'*"e~ were stud- 
ied by considering both relativistic and QCD radiative 
corrections, and predictions were made. With the same 
method of studying the 2S charmonia in the light-front 
quark model, we can study more decay modes, and calcu- 
late the masses of these bottomonia in this paper. From 
an experimental viewpoint, a large amount of bottomonia 
and their excited states could be produced at the forth- 
coming LHC or by the Belle experiment in the near future, 
and they could provide important tests of different predic- 
tions. 

Moreover, heavy quarkonia, especially charmonia and 
bottomonia, act as improtant diagnostic tools to probe 
the properties of the background QCD matter, such as of 
the formation of the quark-giuon plasma (QGP), in heavy 
ion collisions at RHIC and LHC [17] . As has been pointed 
in Ref. [T7], the study of heavy- quarkonium suppression 
at RHIC energy, which might be a signature for the QGP 
formation, calls for the knowledge of the light-front wave 
functions of the quarkonia: /(xx, x^, tq) = ip{x±)ip* {x'^), 
where (p{x±) can be taken as the Fourier transform of our 
light- front momentum space wave functions for the IS or 
2S quarkonia (Eqs. ([7]) and dS])), and f{x±,x'^,To) are the 
essential quantities to calculate the transverse momentum 
distribution of quarks. Thus our wave functions can not 
only help us understand the structure of heavy quarkonia 
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themselves, but also be used as inputs for other physical 
studies. 

This paper is organized as follows. In Sec. [H we de- 
scribe the light-front quark model and the 2S state wave 
function for the quarkonia. In Sec. |31 we present our nu- 
merical results of the decay constants, form factors and 
masses of these charmonia and bottomonia, and compare 
them with experimental data. A brief summary is given 
in Sec. H 



2 Model description 

Heavy quarkonia have been studied by non-relativistic 
treatments [T81IT91I20] . but in some occasions related to 
non-perturbative scales, they have to use model depen- 
dent methods. And as the virtual photon momentum 
increases, the relativistic effects become important. So it 
is useful to study quarkonia in a relativistic treatment. 
Several powerful non-perturbative tools have been devel- 
oped to study the structure and decays of mesons, such as 
the QCD sum-rule technique and the lattice gauge theory. 
The light-front quark model is also an important model to 
do such studies 21,22,23 , and it has a number of salient 
features. Light- front quark model includes some impor- 
tant relativistic effects that are neglected in the traditional 
constituent quark model, and the vacuum in the light-cone 
coordinate is simple because the Fock vacuum is the ex- 
act eigenstate of the full Hamiltoian. Light-front quark 
model has been successfully applied in many investiga- 
tion of hadron structures [SIlEllSllMlESlEaEZllMllMlllSl 

In the light-front quark model, the states of quarkonia 
can be described by the Fock state expansion 



\M) - 



(1) 



and to simplify the problem, we adopt the lowest order of 
the above expansions and take only the quark-antiquark 
valence states of the mesons into consideration. 

The quarkonium wave function in light-front formalism 
is [1I1[22,39 



•0(a;,ki)xi|(x,kx,Ai,A2), (2) 

with the momentum of the struck quark being (xP^ , [to'^-I- 
{xP±+k±)'^]/{xP'^),xP±+k±), and being the helicity 
of the i-th constituent quark. (j){x,'k.±) is the radial wave 
function, and Xm '^■-Lj '^ii -^2) is the light-front spin wave 
function, which is related to the instant-form spin wave 
function by the Melosh-Wigner rotation |40ll411lt^ E5 1 [M l 



xliT) = wdiK + m,)xi (F) - k^xtiF)] 
X\{T) = w,[{kt + m,)x\{F) + kfx]{F)] 



(3) 



where = l/y^2fc+(fcO + m,), fc^^-^ = fc^ ± fc^, fc+ = 
k'^ + k'^ — xjv[ , mi is the mass of the constituent quark, 



and the invariant mass of the composite system = 

^ (k^^ + mD/x + (k_L^ + mD/il - x). The Melosh-Wigner 
rotation is an important ingredient of light-front quark 
model and plays an essential role in explaining the "proton 
spin puzzle" [24 , 25 , 26i.[g71[gl[gl[5ni . In the above formal- 
ism, the Drell-Yan-West (g"*" = 0) frame [^H5] is used be- 
cause only valence contributions are needed in this frame 
when studying the decay of quarkonia. 

For the radial wave function (/>, the harmonic oscillator 
wave function has been adopted to describe the IS state 
mesons [39ll36ll37ll38] , and it can well explain experimen- 
tal data. So we try to go further to use the 2S state har- 
monic oscillator wave function as the approximate wave 
function of the 2S quarkonia. The wave functions of the IS 
and 2S states of the non-relativistic 3-dimensional isotrop- 
ical harmonic oscillator in momentum space are 



1 

■exp( 



(4) 



3.3/4^n)V2 (P' - ^a^;.^)exp(-^),(5) 



where a = fj,uj/h, /i and uj are the mass of the oscillating 
particle and the frequency of the corresponding classical 
oscillator respectively. 

We use the connection between the equal-time wave 
function in the rest frame and the light-front wave function 
suggested by Brodsky- Huang-Lepage PT1[^I55] . for the 
quarkonia with mi = ni2 = niq, 



P 



4:X{1 



(6) 



and we use the prescription in Ref. |46| to extend the non- 
relativistic form wave function into a relativistic one |36) . 
Then we have the corresponding relativistic wave func- 
tions in light-front formalism 



0^^(a:.k.,) = ^^— exp(-— ), (7) 



^ - - 2^ )^^p(-2^^^) 

where /? is the parameter equivalent with a in Eqs. dH) 
and ([S]), and its value can be chosen to fit experimental 
data. The longitudinal momentum k^ = {x — 1/2) -|- 
(rrij ~ m1)/2Ai, one can easily check that this is equiv- 
alent to Eq. ©. The factor ^/dkjd X in the above two 
equations comes from the Jacobian of the transformation 
(a;,k) — ^ (k, fc^), and the normalization factors are from 
the requirement of the normalization of the total wave 
function [47] , 

Using the above formalism and wave functions, we can 
calculate the decay constants and transition form factors 
of the quarkonia [151115] . 

In the V — e+e~ process, the decay constant of the 
vector meson V is defined by 



{'^\3^V{p,S,))=Mvfve^,{S. 



(9) 
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and with the same method as Ref. [SD], we have, for the 
vector quarkonium, 



r- /-dxd^k, 



1 



- x) 



'(a;,k_L) 



2kj_ +mq{M + 2mq) 



(10) 



with the last term being the hyperfine interaction that 
causes the mass splitting between vector and pseudoscalar 
mesons. Here we choose the confining potential (the sec- 
ond term) to be the harmonic oscillator potential rather 
than the linear potential in order to keep consistency with 
our harmonic oscillator wave function for the quarkonium. 
The values of parameters a, b and as were given in Ref |51) 

The mass of the meson is obtained as Mgg — {ip\Hqq\ 
For the 2S quarkonium, we have 



where rriq and is the mass and electric charge of the 
constituent quark of the quarkonium respectively (e^ = 
2/3 for charmonia, and —1/3 for bottomonia). 

In the P — > 77* process, the transition form factor of 
the pseudoscalar meson P is defined by 

{l{p-q)\J''\P{p, A)) = te'Fp^^,, {Q^)e^-P'' p,ep{p~q, X)q,, 

(11) 

and we have the formula for the pseudoscalar quarkonium 



- ^6 



167r3 
x(l — x) 



(12) 



where = — (1 — a:)q_L, and = —q^ = qIi 9 is 
the momentum of the virtual photon. 

The radiative transition form factor between a vector 
meson V and a pseudoscalar meson P is defined by 



{P{p')\J^^\V{p,\)) = iei^y^^p(Q2)£M-p-e^(p,A)p>, 



and we have 

mq{M + 2mq){l - x) + 2(1 - x)kl^ sin^ 6 



(13) 



{M + 2mq)J'k\ + m'^Jk'l + 



,(14) 



where 9 is the angle between k_L and . 

The above quantities are related to the decay width of 
the quarkonium by |48) 



r{V e+e 



3My ' 



r{P ^ 77) = -,ira^Ml\Fp^^^.m\ 



(15) 
(16) 



lv~^lP - ^\i'V^TPW\ 1 7m ) -il') 



2Mv 



We can also calculate the mass of the quarkonium, 
using the QCD-motivated Hamiltonian for mesons |51] 



2^ +k2 



k2 + Vaa, 



'99 ~ V '"g ^ ^ V "'<3 ^ ''qqi (18) 

where k is the momentum of the constituent quark, and 



.9 4q;s 2So - So 9 , , 
Vqq ^a + br^--^ + v' 
3r Smqniq 



(19) 



7b 20asl3 
st^m^ (vector quarkonia). 



(pseudoscalar quarkonia) 



(20) 



and such formula of the mass of the IS quarkonium can 
be found in Ref [52; . 



3 Numerical results 

In our numerical calculation, the parameter /3 in the wave 
function and the mass of the constituent quark rriq were 
chosen to fit experimental data. Since the only difference 
between the vector and pseudoscalar quarkonia that share 
the same energy quantum number (n) is the hyperfine 
interaction term in this model, we choose the same /3 
for them. For charmonia, vcic and /3j/^ {.^^■^^ were fixed 
by Refs. [171152] . and their results are in good agreement 
with experimental data, so we use their values of rric and 
/3j/^ (/3r)c)j s-iid we only have to fix the parameter /3^/ (/3,;^). 
For the bottomoina, we fix all the parameters rrifc, Prif^^-ib) 
and Pr'WnO- 

The parameters of the charmonia are fixed as 

TOc = 1.8GeV, (/3,, J = 0.6998 GeV, 

/?^,(/?^,) = 0.630 GeV, (21) 

and our numerical results of 2S charmonia are listed in Ta- 
ble [TJ The numerical results of IS charmonia can be found 

Table 1. Numerical results (GeV) of 2S charmonia. 



Experiment 16 



Theory 



.e7-(0) 

/^/(V'' e+e 

^< 
F^i- 
F„,. 



>17^7* (0) 
^77* (0) 



0.0392±0.0031 0.0402 

0.1910±0.0057 0.2474 

3.686093±0.000034 3.778 

3.637±0.004 3.637 

< 0.9006 0.6292 

< 0.0590 0.0271 



in Refs. [471152] . The transition form factors F^i ^^i' -y' {^) 
and F,Y (0) are our predictions, and we see from the 
table that they are well below experimental upper limits. 
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We can also obtain the branching ratios of the two decay 
modes using Eqs. ([T6)) and (|T7)) and the total widths of xp' 
and [T^: 

Br(V'' ^ 777c) 3.9012 x 10"'', 
Br(?7^ ^ 27) = 1.0555 x lO"''. (22) 

BES-III collaboration reported very recently the first 
measurement of the branching ratio Br — > 7/7^) = (4.7± 
0.9stat±3.0sys) X 10"^ [53], and our prediction in Eq. ^ 
is in agreement with the preliminary data within error 
bars. The very recent theoretical study in Ref. |8J gives 
777c) = 0.08^^;J53 keV, and converted into the 
branching ratio using the total width of ip' [IH] , it is Br(-0' — >■ 
777^) = (2.7972±1.1) X 10-*. Other theoretical predictions 
for Br(V'' ^ 777^,) fall in a range of (0.1 - 6.2) x lO"" [M]. 

Assuming that 77c and 77^ have equal branching frac- 
tions to KsKtt, Ref. [5S] obtained the experimental data 
^■y-yiVc) = 1-3 i 0.6 keV, using the total width of rj'^ [16 , 
the branching ratio is Br(77^ 27) = (0.9286 ± 0.63) x 
10"''. The theoretical prediction in Ref. pj, converted into 
branching ratio, gives Br(77c —> 2j) = 1.4286 x 10"^. 
We see that our prediction in Eq. ((22|) is close to these 
data. However, the accurate experimental data for this 
decay mode is still not available, and only the upper limit 
Br(r7c — 27) < 5 x 10~* is given \W . Future experimen- 
tal measurements at BES and CLEO may provide tests 
for these predictions of Br('0' — 77;^) and Br(77c ^ ^7). 

For bottomonia, the parameters are fixed as 

mb = 5.1 GeV, PriP^,) = 1-1656 GeV, 
/3y,(/3^,) = 1.1050 GeV, (23) 

and our numerical results of the IS and 2S bottomonia are 
listed in Table [21 Our results give the prediction 1^(77^ — > 

Table 2. Numerical results (in units of GeV) of the IS and 
2S bottomonia. 



Experiment 16 Theory 

fr'{T' e+e") 0.1657±0.0097 0.1944 

fr{T e+e~) 0.2384±0.0044 0.1822 

Fr'^rj^^'iO) -0.0047±0.0009 -0.0047 

Mr' 10.0233±0. 00031 10.0592 

Mr 9.46030±0.00026 9.4087 

M^, 9.3909±0.0028 9.3346 

M^> ? 9.9644 

Fr,''-,^^.^ (0) ? 0.0019 

F^l^^j. (0) ? 0.0020 

-Fr'^^^7-(0) ? -0.1216 

Fr^rj'-y' (0) ? -0.1261 



27) = 0.1494 keV, compared with the prediction given in 
Ref. [T] r(77^ — > 27) = 0.21 keV. Future experiments at 
LHC or by the Belle experiment on 77^ can not only test 
these predictions, but also help us learn more about this 
meson by providing more experimental information about 
it. 



The small values of the experimental data for the branch- 
ing ratio of the mode V{2S) jP{lS) in Table □ and 
Table[2]can be easily understood with our wave functions, 
as the 25' and IS* wave functions are orthogonal to each 
other and their overlap in Eq. (jl3p is suppressed. 

Although nonrelativistic models have provided efficient 
and powerful theoretical tools to handle various problems 
related to hadron structure, the relativistic models have 
been successful in many investigation of hadron struc- 
tures |3ii[n3[S3[Mi[33[3a[szi[Mii2i[23[ia[2zi[ia[2a[so]. it 

is thus necessary to make an estimate of the effect due to 
nonrelativistic to relativistic treatments [55]. For simplic- 
ity, we assess the non-relativistic to relativistic effects by 

letting k^^ of y^k^ -I- in the expressions of the decay 

constants and form factors to be zero. For examples, af- 
ter this procedure, we have F^'^^^^.(O) = 0.1167 GeV, 
compared with F^/^^^^.(0) = 0.0402 GeV/0.0392 GeV 
from the relativistic treatment / experimental data, and 
fr'OT' ^ e+e") = 0.2159 GeV compared with fr'{r' 
e+e-) = 0.1944 GeV/0.1657 GeV from the relativistic 
treatment / experimental data. We see that the relativis- 
tic treatment is needed to describe the experimental data 
well. 



4 summary 

In this work, we studied the 2S quarkonia ip' , 77^, T' and 7;^ 
in light- front quark model. Similar with the IS harmonic 
oscillator wave function that was commonly used as the 
IS quarkonium wave function in light- front quark model 
studies, we tried to use the 2S harmonic oscillator wave 
function as the 2S quarkonium wave function. The decay 
constants and transition form factors of these quarkonia 
are calculated. Using the QCD-motivated Hamiltonian for 
mesons, we also calculated masses of these quarkonia. Our 
numerical results of these quantities are in agreement with 
experimental data. Predictions of transition form factors 
and masses of these quarkonia are made, and these predic- 
tions can be tested by future experiments. The IS and 2S 
wave functions could also be used as inputs to study other 
problems such as the "p7r puzzle" and the suppression of 
heavy-quarkonia at RHIC energy. 
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ral Science Foundation of China (Grants No. 11021092, 
No. 10975003, No. 11035003, and No. 11120101004) and 
by the Research Fund for the Doctoral Program of Higher 
Education (China). 
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